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Abstract. We show that a simple quantum-mechanical model, put forward by Stone 
sometime ago, affords description of site magnetoelectricity, a phenomenon which takes 
place in crystals (and molecular systems) when space inversion is locally broken and 
coexistence of electric and magnetic moments is permitted by the site point group. 
We demonstrate this by identifying a local order parameter, which is odd under both 
space inversion and time reversal. This order parameter (a magnetic quadrupole) 
characterises Stone's ground state. Our results indicate that the model, extended to a 
lattice of sites, could be relevant to the study of electronic properties of transition-metal 
oxides. A generalisation of Stone's hamiltonian to cover cases of different symmetry is 
also discussed. 
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In 1986 Stone reported a study of the Hamiltonian 

H = jjL^~fin-(T, (1) 

which he performed by path-integral and by conventional quantum-mechanical 
techniques (Stone 1986). The conventional method considers the basis set of spinor 
spherical harmonics 

as suggested by [iJ, J]_ = 0, with J = L + jcr; L stands for the angular momentum 
operator that generates rotations of n = r/r; cr are Pauli matrices and I is the 
rotor moment of inertia; [, ]_denotes a commutator and Cji^i;j2,m2 ^ Clebsch-Gordan 
coefficient. On account of the property (Stone 1986, Varshalovich et al 1988) 

n ■ (t\j ± i, jm) = -|j T l,jm) , (3) 

the Hamiltonian H is readily diagonalised in this basis. In the large-/i limit, the ground 
state is found to be 

1^)- = 72 i\3 + l'-^'"^) - 1^' - 1'^'"^)) ' (4) 
with energy E = —fi. 

Stone's motivation was to provide a simple quantum-mechanical example in which 
the Berry phase gives rise to Wess-Zumino terms. Indeed, for large /i, Eq. ((T)) describes 
the motion of a constrained spin, which is equivalent to motion of a charged particle 
about a magnetic monopole (Leinaas 1978). 

Equation (^J) has been interpreted as a model for a solenoid, which is rotating about 
its centre of mass where a spin-| particle is placed. When fi is small, the solenoid and 
the particle would spin independently. As ^ becomes large, the spin will become slaved 
to the direction of the solenoid (Stone 1986, Aitchinson 1987). Notice that this physical 
picture leads to the coupling B ■ cr, which is space and time even. {B stands for the 
magnetic field generated by the solenoid.) Such a symmetry should be contrasted with 
that of n ■ cr, which is space and time odd. 

The current Letter will demonstrate that the model defined by Eq. affords 
description of a different effect: site magnetoelectricity. Such a phenomenon occurs in 
crystalline and molecular systems when space-inversion symmetry is locally broken and 
co-existence of electric and magnetic moments is permitted by the pertinent site point 
group. An effective magnetoelectric interaction between these two moments would be 
described by Eq. provided we identify n with a unitary electric-dipole moment. 
(The electric charge e is merged into /x.) This new interpretation of the model does not 
affect its dynamical regimes (Stone 1986), which remain those of the rotating solenoid 
with B replaced by n. 

As shown by Goulon and his collaborators (Goulon et al 2000, 2002), microscopic 
magnetoelectric behaviour of crystals can be investigated using near-edge absorption of 
X rays, which implies excitations of inner-shell electrons to empty valence states. As is 
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known, this experimental technique is site selective, a feature resulting from the tuning 
of x-ray energy at a given inner-shell threshold. Sensitivity to the long-range order of 
local magnetoelectric order parameters is obtained by recording dichroic signals which 
stem from an interference between electric-dipole and electric-quadrupole transitions. 
As a consequence, scalars (e.g. n-cr) are not probed by these experiments, which detect 
the long-range order of local (on-site) magnetoelectric order parameters represented by 
one-particle irreducible tensors of rank 1,2 and 3. One set of these order parameters 
specifically serve our purposes: the magnetic quadrupoles (rank-2 tensors) \. In the 
LS-coupling scheme, they read 

M^^^^[n,Lf\ MP^[n,Sf\ (5) 

as shown by recent theoretical work on x-ray dichroism and resonant scattering in 
noncentrosymmetric crystals (Carra et al 2003, Marri and Carra 2004). (The symbol 
[,]'^'^) denotes Clebsch-Gordan coupling of irreducible tensors; the spin operator is 
defined hj S = |cr; Q^^) = [L^L^^^ denotes an orbital quadrupole.) Inspection of 
magnetic nature identifies orbital, spin and spin-orbital elements in the set. All tensors 
in (0) are space and time odd, being thereby invariant under the combined action of 
these transformations. 

A set of vector order parameters will also be considered in connection with Stone's 
model. Its elements are defined by (Marri and Carra, 2004) 

n, Ps^^lxS and = -^-^[[n, Lf\ Sf\ (6) 

in LS coupling. These irreducible tensors have polar (electric) symmetry, i.e., they are 
space odd and time even §. 

The order parameters (jSHHl) are defined in the second quantisation formalism. For 
example, 

and similarly for the others. Here, c|^^ and c^^^ denote fermionic creation and 
annihilation operators, respectively. 

As stated, the main purpose of the current work is to discuss magnetoelectric 
properties of Stone's Hamiltonian in the large n limit. To this end, we will show that 
the symmetry property displayed by the scalar n-cr, when acting on the spinor spherical 
harmonics |j±|, jm) [Eq. Q], extends to the irreducible tensors defined by expressions 

I A simple mechanical model of a magnetic quadrupole is provided by two parallel coils run through 
by opposite currents. 

§ The LS tensors (|5I6|I are given in compact forms. In this representation irra- 
tional prefactors appear. They can be removed by recoupling transformations. We 
find [5, [i,n](2)](i) = (^s ■ Ln + L S ■ n) , and similarly for [5, [^l, i]^^^]^^^; further- 

more, [S,-z[f2i,i](2)](2) ^([SxJli,i](2) + i[S,fliXi](2)) and [5, [n, Q(2)](3)](2) ^ 

{3[S, [n, Q(2)](2)](2) + 6[[5, Q(2)](2)^„](2) _ 3[_5 . ^ g(2)](2) „ 5^^^ [g(2) ^ (2) ] (2) ) ^ 



{a'\{l'm'\{MPUlm)\a)cl 



+ h.c. 



(7) 
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(iniini), leading to relations whereby the magnetoelectric behaviour of \g)- is readily 
inferred. 

The basis set \j ± ^,jm) provides a convenient framework for describing parity- 
breaking electron hybridisation (e.g. pd mixing in transition- metal oxides), in the jj 
coupling scheme ||. It is thus clear that to generalise Eq. (jH)) to higher-rank order 
parameters we must first determine the form of the corresponding tensors in jj coupling. 
This is accomplished by resorting to the theory of LS— >jj transformations (Edmonds, 
1974). In the case of space-odd irreducible tensors, such transformations contain matrix 
elements between the states \jm) and \ j'm'), with j = / ± ^ and j' = /' ± |. Solving the 
corresponding equations for j' = j, as demanded by [H, J]„ = 0, provides the required 
jj-coupled order parameters, which will appear as linear combinations of space-odd LS- 
coupled irreducibe tensors. (Notice that n-cr/2 = n- J, as n- L = 0; no transformation 
is needed in this case.) Technically, the derivation is laborious as it implements angular- 
momentum recoupling methods (Racah calculus). Thus, for convenience of the reader, 
we will first state our results, then illustrate their physical content and, at the end, 
discuss mathematical aspects of the formulation. 

The required jj-coupled magnetic quadrupole is found to be 

Mfil'J) = I (^) (^) MPil',l) + IM?{1',1) 

+ lMPil'J)-lM^Pil',l) (8) 

yielding, after some algebra, 

E Mf{l'J).\j±l,jm) = -'-^^!^^mi\j^l,jrn), (9) 
i,i'=i±i 

with z the quantisation axis. This result generalises Eq. (jS)) and shows that \g)_ is an 
eigenstate of the jj-coupled magnetic quadrupole operator; in other words, the large- 
/i ground state of Stone's hamiltonian is magnetoelectric. Our conclusion is further 
supported by what follows. Consider the jj-coupled operator 

Pj{l', I) = n{l', I) + Ps{l', I) - 2Pt{1', I) . (10) 

We have 

J2 PAl',l)o\j±hjm) = -m\jThjm), (11) 
i,i'=i±i 

showing that \g)^ is an eigenstate of the jj-coupled unitary electric-dipole moment, with 
eigenvalue m. \g)_ is thus characterised by the simultaneous presence of an electric and 
a magnetic moment in a parallel (as expected) configuration. [Reversing the sign of 
the coupling constant in Stone's model (/i — /i, large /i) would change the ground 
state to \g)+ = (|j + \,jm) + \ j — ^,jm)), which is characterised by an antiparallel 
alignement of the moments and by a magnetic quadrupole with opposite sign.] 

II As shown below, in this basis, parity breaking hybribisation displays rotational symmetry [SU{n), 
in the general case]. Rather than by ordinary \jm), irreducible representations are spanned by 
Tjdi + 5, j?™) ± \j - for ranks and 2, and by ^(e^/^lj + 5, jm) ± e-*''/'*|j - ^Jm)) 

for rank 1. 
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Thus, we are led to the conclusion that, in the strong coupling limit, the 
magnetoelectric interaction can be viewed as a problem of a magnetic moment 
constrained by an electric field. In turn, this is equivalent to the problem of a charged 
particle moving in the field of a magnetic monopole, both classically and quantum 
mechanically (Leinaas 1978). This seems to tally with recent work on magnetic 
monopoles in crystal momentum space (Fang et al 2003). 

According to our findings. Stone's model provides a good starting point in the 
study of interactions between (local) electric and magnetic moments in crystals. For 
this purpose, an extension of the model to a lattice of sites is now needed. Such a 
model, characterised by an order parameter which violates space inversion and time 
reversal, could be relevant in the analysis of electronic properties of transition-metal 
oxides. 

In certain site point groups (Cracknell 1975), magnetoelectric interactions result in 
configurations where electric and magnetic moments are mutually perpendicular. Such 
interactions arise from toroidal distributions of currents and are described by anapole 
moments. In LS coupling, the set of anapolar order parameters reads (Carra et al 2003, 
Marri and Carra 2004) 

VLl = ^{n X L — L X n) , ^Is = n x S , 

(12) 

The required jj-coupled anapole takes the form 

Oj(r, I) = ^f^L(r, /)-(/ + /' + i)ns{i\ i) + j^.^rii', i) , (is) 

giving 

(Oj),|j±i,jm)= ^j{l\l)z\3±\,3m) = Ttm\jT\,3m). (14) 

l,V=l±l 

\g)- = -^{e^'^^'^lj + \ —e~'^^/^\j — \ ,jm)) is therefore an eigenstate of the jj-coupled 
anapole operator. It is readily shown that |^)_ is the large-// limit ground state of the 
hamiltonian obtained by replacing n ■ cr with z ■ fij in Eq. (Q). Notice that \g)- is 
also an eigenstate of a jj-coupled pseudodeviator, a polar rank-2 tensor, with eigenvalue 
(j + l)]/V6. This new hamiltonian thus displays x-ray natural circular dichroism. 
[The derivation of this result will not be given for lack of space. LS pseudodeviators 
have been discussed by Marri and Carra (2004).] 

We conclude with some remarks concerning technical aspects of our derivation. 
Generalising the concept of coupled double tensor (Judd, 1967), we define 

W^f'^'{l'J)= Yl ^:Lq:,|.^^'i';a4AV'£a. + h-C. (15) 

and 

Vp>il', l)=Y Cf^,,Jl,rr^^,m + h.C. , (16) 
m,m' 

where c^^^ = (— 1)' '^^^ '^c^-a-o- ^ijm ~ (~l)''~'^'^z j-m' ^^^^ creation and 
annihilation operators transform as the components of irreducible tensors. [Notice 
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that w^^^^^ {I' , I) are spinless, whereas w^^^^^{l' ,1) depend on spin.] The importance 
of w^^y^^{l', I) lies in the fact that all one-electron LS order parameters defined by dSJEJ 
can be expressed as multiples of them (Wigner-Eckart theorem). For example, 

In a similar way, one-electron jj order parameters can be expressed as multiples of 
(/',/)• The tensors w^^^^^{l\l) and v^^'^^^{l',l) are related by a standard LS— >jj 



V 



transformation (Edmonds, 1974), which reads 



I I' X 

i 1 y }v^^'^^^{l'J), (18) 
j j' z 



with [a, b] = (2a = l)---(26+l). In the case of magnetic-quadrupole order parameters, 
Eq. (|T8|l yields a system of four equations. By solving this system for f = j and /' = /±1, 
we find 

Mf{l\l) =r.h.sof Eq. © 

= -1(2/ + l)(2r + 1) {[n, Jf'+h^-^'5v,-^ + [n, J](''-i'+^)^ . 

Vector order parameters lead to systems of three equations. Their solutions for 
j' = j and /' = Z ± 1 read 



and 



3(/+r+l) 



n , 



4(2;+l)(2r+l) 

i+i'+i 



' 2 ''^2 A 



(19) 
(20) 



providing a full definition of Eqs. ^ and In Eq. (j2I|, 

^^?'(^'.0 = 5^(r,jW|(Vn- Vj,) X J + Jx (Vn- Vt,)|/,jm) 



where = —in x L ^. 

In the case of anapolar order parameters, LS Judd's tensors are defined as 

and similarly in jj coupling. 



(21) 



(22) 



I am grateful to M. Fabrizio, A. Jerez, E. Katz and T. Ziman for discussions and a 
critical reading of the manuscript. 

^ A natural choice for the anapole in jj coupling would be (n x J — J x n)/2 (Dothan and Ne'eman, 
1966). However, this operator cannot be employed in our case as its matrix element vanishes for j' — j. 
The vector appearing in Eq. removes this drawback (Jerez, 2003). 
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